full_join(
left_join(EB_preference, EB_rrs) %>%
mutate(reporting_rate = paste0(substr(
formatC(reporting_rate * 100, format = "f", digits = 4), 1, 5), "%")) %>%
mutate(EB_preference_est = if_else(
is.na(EB_preference_est), "NA", if_else(
is.infinite(EB_preference_est), "Inf", substr(
formatC(EB_preference_est, format = "f", digits = 4), 1, 5)))) %>%
mutate(EB_preference_lci = if_else(
is.na(EB_preference_lci), "NA", if_else(
is.infinite(EB_preference_lci), "Inf", substr(
formatC(EB_preference_lci, format = "f", digits = 4), 1, 5)))) %>%
mutate(EB_preference_uci = if_else(
is.na(EB_preference_uci), "NA", if_else(
is.infinite(EB_preference_uci), "Inf", substr(
formatC(EB_preference_uci, format = "f", digits = 4), 1, 5)))) %>%
mutate(EB_preference_ci = if_else(
EB_preference_lci == "NA" | EB_preference_uci == "NA", "NA",
paste0("[", EB_preference_lci, "; ", EB_preference_uci, "]"))) %>%
mutate(EB_preference = if_else(
CI_calculable == TRUE & EST_calculable == TRUE,
paste(EB_preference_est, EB_preference_ci),
"No estimate")) %>%
select(scientific_name, "EB_reporting_rate" = reporting_rate, EB_preference),
left_join(BT_preference, BT_rrs) %>%
mutate(reporting_rate = paste0(substr(
formatC(reporting_rate * 100, format = "f", digits = 4), 1, 5), "%")) %>%
mutate(BT_preference_est = if_else(
is.na(BT_preference_est), "NA", if_else(
is.infinite(BT_preference_est), "Inf", substr(
formatC(BT_preference_est, format = "f", digits = 4), 1, 5)))) %>%
mutate(BT_preference_lci = if_else(
is.na(BT_preference_lci), "NA", if_else(
is.infinite(BT_preference_lci), "Inf", substr(
formatC(BT_preference_lci, format = "f", digits = 4), 1, 5)))) %>%
mutate(BT_preference_uci = if_else(
is.na(BT_preference_uci), "NA", if_else(
is.infinite(BT_preference_uci), "Inf", substr(
formatC(BT_preference_uci, format = "f", digits = 4), 1, 5)))) %>%
mutate(BT_preference_ci = if_else(
BT_preference_lci == "NA" | BT_preference_uci == "NA", "NA",
paste0("[", BT_preference_lci, "; ", BT_preference_uci, "]"))) %>%
mutate(BT_preference = if_else(
CI_calculable == TRUE & EST_calculable == TRUE,
paste(BT_preference_est, BT_preference_ci),
"No estimate")) %>%
select(scientific_name, "BT_reporting_rate" = reporting_rate, BT_preference) %>%
left_join(., LUT_taxonomy %>%
select(scientific_name = BT_scientific_name, EB_replacement_name)) %>%
mutate(scientific_name = if_else(is.na(
EB_replacement_name), scientific_name, EB_replacement_name)) %>%
select(-EB_replacement_name)
) %>%
arrange(scientific_name) %>%
replace(is.na(.), "Not in dataset") %>%
left_join(., ebird_taxonomy) %>%
select("Scientific Name" = scientific_name,
"Common Name" = common_name,
"eBird Reporting Rate" = EB_reporting_rate,
"eBird Preference Ratio" = EB_preference,
"BirdTrack Reporting Rate" = BT_reporting_rate,
"BirdTrack Preference Ratio" = "BT_preference") %>%
write_csv("SupplementaryTable1.csv")